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Abstract 



C^^ I Two locally positive expressions for the gravitational Hamiltonian, one using 

f^ ^ 4-spinors the other special orthonormal frames, are reviewed. A new quadratic 

(N 

(^ ' 3-spinor-curvature identity is used to obtain another positive expression for 

o' 

^ , the Hamiltonian and thereby a locaUzation of gravitational energy and pos- 

Qv I itive energy proof. These new results provide a link between the other two 
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methods. Localization and prospects for quasi-localization are discussed. 
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I. INTRODUCTION 

For asymptotically flat gravitating systems total energy is well defined and must be non- 
negative. Each new positive total energy proof (e.g. [jl|) offers some more insights. Con- 
cerning the localization of the total energy, although the equivalence principle forbids a true 
local gravitational energy density yet a suitable "quasi-localization" is desirable . A good 
candidate for a gravitational energy density is the Hamiltonian density. For asymptotically 
fiat Einstein gravity the Hamiltonian density has the general form 

H{N) = j d^x2N^'G\+ Ib 

Sx Nn + N^nk+ iB , (1.1) 



which includes a boundary term at spatial infinity. On a solution the spatial integral van- 
ishes, the value of the Hamiltonian, — IGvrGA^'^p^, comes from the integral of the boundary 
term over the 2-sphere at spatial infinity and determines the total 4-energy momentum p^. 
The integrand of the boundary term, S, is only well defined up to 0{r~'^), moreover we have 
the freedom to choose the lapse N and shift N^ . Together these allow a certain latitude 
which can be exploited to obtain a locally non-negative Hamiltonian density. Indeed, such 
a form can be achieved in more than one way. 

II. THE 4-COVARIANT QUADRATIC SPINOR HAMILTONIAN 

The first constructions of this type [Q were done in the wake of the Witten positive 
energy proof. It was shown that the Hamiltonian density for Einstein gravity could be 
expressed as a 4-covariant quadratic spinor 3-form: 

H(V') := 2{Z}(^757) AD^-D^A D{^,^i/j)} . (2.1) 

This remarkable result follows from: (i) the identity 

H(^) = 2mG''^ri, + d{V^757 AD^Ij + /^(V^TsT)^ 

-V^D(757V') +^?A (757^)} , (2.2) 



where iV^ = i/jj^^tp (for conventions see the appendix) reveahng that ?i(V') contains the 
appropriate projected components of the Einstein tensor (needed to generate the equations 
of motion) up to an exact differential (which does not change the variational derivatives); 
and (ii) the observation that 6 JTilip) has an asymptotically vanishing boundary integral 
since 7i{4') is asymptotically of order 0{r~*). 

The Hamiltonian density ( |2.1|) can be decomposed, with respect to the normal to any 
spacelike hypersurface, into positive and negative definite parts: 

n{^) ~ Aig'^'DaiJ^D,^- I ^Da^ n r/o , (2.3) 

and is locally non-negative if ip satisfies the Witten equation (or certain modifications 
thereof) 

^Dai^ = , (2.4) 

thereby permitting a non-negative "localization" of gravitational energy. 

This mathematically elegant form of the gravitational Hamiltonian has several virtues, 
in particular it (i) is manifestly 4-covariant, (ii) shows that total 4-momentum is future 
timelike, (iii) can be evaluated on a spacelike surface extending to future null infinity thereby 
showing that the Bondi 4-momentum also is future timelike. However it also has some 
liabilities, in particular (a) the spinor field is physically mysterious, (b) there is no direct 
relation to the customary variables, (c) it yields an unintuitive energy localization. (For 
the Schwarzschild solution in isotropic cartesian frames ip = {1 + m/2r)~'^ip^ j^ solves the 
Witten equation; using this result in the boundary integral yields 1/8 of the total mass- 
energy inside the horizon and 7/8 outside.) Consequently other Hamiltonian based positivity 
proofs/localizations were sought and found. 

III. THE SPECIAL ORTHONORMAL FRAME APPROACH 

Another approach H| used orthonormal frames and exploited their rotational gauge free- 



dom. The ADM Hamiltonian (|1 . 1|) in an asymptotically cartesian frame has the form: 



H{N) = fd^x Ar(^-5(7r'""vr„„ - ^n^)-giR) 

+27i^k'^mN' + idSk iV5,\^r«^ . (3.1) 



We choose N'' = 0, use the divergence theorem to ehminate the boundary term, parameterize 
the metric with orthonormal frames, spht the connection coefficients algebraically into a 
symmetric tensor qat, a vector qc := —T'^ca and a scalar q := e'^'^'^Tabc, and use the special 
orthonormal frame (SOF) ^ rotational gauge conditions: 

qk = 4(9fc In $, q = constant, (3.2) 

to obtain the Einstein Hamiltonian (i.e., energy) density in the form 

H{N) = 8gh^'^dn{N<^-')d^^ 

+iV{^-^(7r»V„5 - i^r^) + g^q'^'qab - h')} • (3.3) 

This expression is good for both compact spatial surfaces (in which case q is a non-vanishing 
constant) and for asymptotically flat spatial surfaces (in which case q vanishes). For the 
latter case total energy is well defined; a suitable choice of the lapse gives a positive total 
energy proof. 

Many choices for the lapse give a positive local energy density, in particular A^ = $", (a > 
1). An especially attractive choice is A^ = $ which leads to the gravitational energy density 

n{^) = '^{g-^TT'^'n^b - y)+9Hq''\b - ig')} , (3.4) 

and the value 

E = (levrG)-^ / ni<l>) d^x = {2nG)-^ I g^V'^^dSk , (3.5) 

for the amount of energy localized within a volume V bounded by a surface S. The value 
is non-negative for asymptotically fiat maximal spatial hyper surf aces. The gravitational 
potential $ (generalized Newtonian potential) satisfies the generalized Poisson equation: 

8^5 A$ = n{^) + 167r<!?g^Gp . (3.6) 



(This is just the Hamiltonian constraint and is essentially the scale equation of the usual 
conformal approach to the Einstein initial value constraints see, e.g., Choquet-Bruhat and 
York § ). 

This SOF Hamiltonian has certain virtues, especially (i) the gauge conditions are con- 
formally invariant so SOFs are closely related to the usual variables of the standard initial 
value constraints, (ii) the oscillating physical modes are apparent in the SOF gravitational 
energy density, and (iii) the energy localization is physically reasonable; in particular, all of 
the mass of the Schwarzschild solution (note: $ = (1 + r7i/2r)^^ for the isotropic cartesian 
frame) is within the horizon. (Moreover there is some freedom here; the choice A^ = <l>^ 
produces the same 1/8 inside the horizon as the 4-covariant spinor Hamiltonian for the 
Schwarzschild solution.) However the SOF approach also has certain liabilities, in particu- 
lar (a) the expression concerns only energy, it gives no restraint on the momentum, (b) the 
energy is guaranteed to be locally non-negative only for maximal spacelike hypersurfaces, (c) 
a maximal spacelike hypersurface cannot be extended to future null infinity so this approach 
cannot give the Bondi mass-energy. 

IV. A NEW 3-SPINOR PROOF AND LOCALIZATION 

A link between the special orthonormal frame approach and the 4-covariant quadratic 
spinor form of the Hamiltonian has now been found in terms of a new Hamiltonian based 
gravitational energy positivity proof and localization which uses 3 dimensional spinors. 

The key is a new spinor identity (see appendix) 

2[V{fUa) AVip- Vf'' A V{icrip)] 

= dB- (y^V)^"' A U , (4.1) 

where 

B := (fHa A V^p — Lp^V{iaLp) 

+V{ipHa)ip + (Vv?"^) A iaip . (4.2) 



Using this identity we replace the scalar curvature term, NRg^^'^cPx = Nfl"-^ A (at, and the 
boundary term in the ADM Hamiltonian (|3.1| ) with the left hand side of ([4.1| ) . The Einstein 
Hamiltonian (with A^ = v^V; ^'' = 0) can then be written as 
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if(Ar) = yrf3^(y,V)^-i(vr'"X 

+2[V(v5"^ia) A v? - Vv9^ A "^{icyip)] . (4.3) 



An important property of the 3-spinor Hamiltonian ( |4.3| ) is that no additional boundary 
term at infinity is needed. As Regge and Teitelboim |Q] have nicely explained it is necessary 
that the boundary terms in the variation of the Hamiltonian vanish asymptotically. To 
verify this property for ( |4.3|) we need only check the variation of the quadratic spinor terms 
(see appendix): 

(5[V((^^^a) A V(^ - V(^^ A W{ia^)] 

+d{6{(pUa) A Vv? - 6(p''\/{ia(p) 

+V{ipUa)6f + (Vv5^) A 6{iaip)}. (4.4) 

Since cp ~ constant + 0(l/r) and Sip ~ 0(l/r) the boundary term falls off as 0{r~^), 
therefore vanishing at spatial infinity for aymptotically fiat initial data. 

On a maximal hypersurface vr = 0, the kinetic terms in the Hamiltonian ( [4.3| ) are non- 
negative. The quadratic spinor terms can also be made non-negative. Since the torsion 
vanishes the spinor terms reduce to 

2[V{(pUa) AVip- Vip^ A V{taip)] 

= iig'^'Va^^Vk^ - V,<^VV^Vb<^]C. (4.5) 

Hence the spinor terms are non-negative for any asymptotically constant ip satisfying the 
3-dimensional Dirac equation 



(j'^Va^ = . (4.6) 

This is a linear elliptic equation similar to the Witten equation; essentially the same argu- 
ments show that unique solutions exist. Hence the local energy density is non-negative on 
asymptotically spacelike maximal slices. 

This 3-spinor Hamiltonian approach by itself has most of the liabilities of the other two 
approaches: (a) it yields the same sort of unintuitive energy localization as the 4-covariant 
spinor expression (indeed they have identical values for the Schwarzschild solution but differ 
when Kab 7^ 7^ tp'-f'^ip) , (b) the expression concerns only energy, it gives no restraint on 
the momentum; (c) the energy is guaranteed to be locally non-negative only for maximal 
spacelike hypersurfaces, (d) the maximal spacelike hypersurface cannot be extended to future 
null infinity so it cannot give the Bondi mass-energy. However, for the 3-spinor Hamiltonian 
some of the other 4-spinor liabilities are not so severe since in this case (i) the spinor field is 
not so mysterious, for (ii) there is a relation to the customary variables via the SOF variables 
as we shall show below. Indeed the principal virtue of this approach is that it relates the 
other two methods we have discussed. 

V. RELATIONS BETWEEN THE METHODS 



The 3-spinor Hamiltonian expression (|4.3| ) is intermediate between the 4-covariant spinor 



Hamiltonian ( |2.1|) discussed above and the SOF Hamiltonian (|3.3| ) . 

On the one hand it can be extracted as a piece of the 3+1 decomposition of the 4- 
covariant spinor Hamiltonian. The orthonormal frame components of the metric compatible 
4-connection project into the components of the 3-connection and the extrinsic curvature 



K^. = -r°". hence 



= Ve^ + lK\^iolb]^. (5.1) 



Consequently the quadratic in Dip Hamiltonian decomposes into: (a) quadratic terms in K""}, 
along with quadratic terms in V^ which are essentially the 3-spinor Hamiltonian density 



( [4 .31 ), (b) linear terms in K"-^, (they are of the form 27r™'cVmA^'^ where N'^ = ip^^^^lj) which 
represent the momentum constraint. (Note that the 3-spinor method decouples the spinor 
field from A^'^; this has both advantages and disadvantages.) 

On the other hand the 3-spinor Hamiltonian not only resembles the SOF approach in (i) 
using a vanishing shift, (ii) considering the kinetic terms separately, (iii) relying on maximal 
slices, and (iv) replacing the potential terms by an expression using different variables, but, 
moreover, there is a close relation between the SOF variables and spinor fields via solutions 
to the 3 dimensional Dirac equation ( [4.6|) . 

Indeed the 3-dimensional Dirac equation explicitly depends only on the parts of the 
connection which appear in the gauge conditions ( |3.2|) : 

= a'^ip^c - \qb(y^^ + \iq^ ■ (5.2) 

An asymptotically constant solution to a°''Va<p = can be factored into a magnitude and a 
unitary transformation which determines an SOF [Q . Conversely, expressed in terms of an 
SOF the Dirac equation reduces to (j'^da^^^f = 0, hence (f = $^V^const- 

VI. LOCALIZATION AND QUASI-LOCALIZATION 

Our considerations have been concerned with obtaining a positive localization of the 
total energy by finding a good expression for the Hamiltonian density. Each localization 
depends on the solution to an elliptic equation, which, in turn, depends on the values on 
the boundary of the region. Since the boundary is at spatial infinity we can simply choose 
suitable constant values as the physically appropriate boundary conditions. 

Beyond distributing the total gravitational energy, there is considerable interest in "quasi- 
localization" , i.e., determining the amount of energy in a finite region without reference to 
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what is outside. The expressions we have discussed could also be used for a finite region. 
Then each of the (locally positive) Hamiltonian densities provides a quasi-local energy. The 
value of the positive quasi-local energy can be obtained from the associated boundary in- 
tegral. The quasi-localization will depend on the choice of boundary values on the finite 
2-surface bounding the region. We, however, do not know how to decide which values on 
a finite boundary are a physically good choice. Several "quasi-localizations" based upon 
4-covariant spinor expressions like (|2.1|) and ( p.2|) have been investigated by others [llO[ . 



Their methods of choosing boundary values for the 4-spinor field could also be adapted to 
our orthonormal frame or 3-spinor fields. Canonical investigations associated with a finite 
region, with particular attention to the possible boundary terms and their relation to what is 



held fixed on the boundary, have been done for the standard variables [|TT|. Such a study of 
the spinor or SOF parameterized Hamiltonian should provide some guidance for the choice 
of appropriate boundary values for finite regions. 

Nontrivial examples of the localizations produced by the 4-spinor, SOF and 3-spinor 
techniques, e.g., for the Kerr solution, would be instructive. However, as noted, the local- 
izations depend on solving an elliptic system of equations, essentially the Dirac equation. 
Unfortunately, aside from the aforementioned spherically symmetric case, there are hardly 



any known exact solutions for the Dirac equation in curved spacetime ||12 . 

Forgoing direct comparisons for actual solutions we can compare the expressions by 
considering desirable properties. We know of no gravitational energy localization method 



which is satisfactory. One list [|T3| for example, requires (i) zero for fiat spacetime, (ii) the 
standard value for spherical solutions, (iii) the ADM value for an asymptotically fiat slice, 
(iv) the Bondi value for an asymptotically null slice, (v) the irreducible mass for the apparent 
horizon, (vi) positive and monotonic. Of the methods considered here, the quadratic 4-spinor 
expression certainly fails (ii) & (v). The SOF Hamiltonian satisfies the positivity requirement 
(vi) only on maximal slices, while the maximal slice restriction precludes satisfying (iv). The 
new 3 spinor technique has all of these failings. 



VII. CONCLUSION 

We have presented a new positive total energy proof for asymptotically flat Einstein 
gravity. The proof uses a 3 dimensional spinor parameterization of the Hamiltonian and a 
new 3-spinor-curvature identity. What insights has it yielded so far? 

As a proof, considered on its own, this method has no advantages and indeed is less 
general than some other known proofs. More interesting is the fact that it provides another 
independent method for obtaining a positive localization of gravitational energy; yet again, 
as a localization method, it has no apparent advantages. 

Probably the most interesting thing is that it provides a link between two other Hamil- 
tonian based proofs and their associated localizations. At the very least this link connects 
the somewhat mysterious Witten spinor field proof and localization to the more usual type 
variables. 

Perhaps this link will play an essential role in finding a modification of our expressions 
into a better Hamiltonian density — one which permits a positive energy proof and gravita- 
tional energy localization combining the virtues of the 4-covariant spinor Hamiltonian and 
special orthonormal frame approaches without the liabilities. Such an expression, compli- 
mented by a good choice of boundary values for finite regions would provide a physically 
reasonable quasi-localization of gravitational energy. 
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APPENDIX: CONVENTIONS AND IDENTITIES 

Our 4-dimensional conventions are: metric signature (— 1, +1, +1, +1), orthonormal 
coframe 6"^, unit volume element rj, unit 3-form basis ?7^ = *^^. Spinor conventions 
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75 = 7°7W^7^ 

It^lv + Ivlt, = -^g^u = 2diag(+l, -1, -1, -1), 

7 = 7^6*^ and Dip = dip — ia;^'^7[^7j,]-?/' is the covariant differential. 
Tlie 3-diniensional spinor conventions used are: a = acO'^ wliere 

witli (Tab = lic^a, o-fe], SO (TcO-ab + (Tahcrc = Kabc, and (ab = CabcO"" , wliere (abc is tlie 3-dimensional 
Levi-Civita tensor witli C123 = +1- 

Tlie identity connecting the 3-diniensional scalar curvature to the spinor expression in 
the Hamiltonian can be verified as follows: 

2\SI{if^i(r) A Vv? - Vv?^ A V{i(rif)] 

= dB - [-(/^■^io-VV + V^(v9'^io-)v9 - if^V'^{i(yif) + {V'^if^)i(yif\ 
= dB - if]"^ A [-v^^io-o-afev? - i^^iaaabi^ - ^p^ia-abO-^p - ip'^iaabO-ip] 
= dB + fi"^ A [v^iaaab + i(rab(r)^] = dB - (v^V)^"' A Cab, (Al) 

where B := (^''"io- A Vip — ip'^Viiaif) + V(v9''"i(T)v9 + (Vy?^) A ■io'(/9. 

Similarly, we calculate the variation of the quadratic spinor terms 

5\S/{if'^i(j) AVip- Vip^ A V{icrip)] 
= [V6{ipHa) A Vv? - V^v?"^ A V{i(rip) + V{ipU(r) A V^v? - (Vv?"^) A V5(io-v9)] 

+ [(5V)(v9^m) A Vv? - {6V)(p^ A V(mv9) + V(v9"^io-) A (5V)v9 - Vv?^ A {6V){i(rip)] 
= d[6{(pUa) A Vv? - (5v9^V(mv9) + V{ipUa)5ip + (Vv?"^) A 6{ia(p)] 
-[-6{ipUa) A VV - Sip^V'^{icrip) + V^{ipU(r)6ip + VV^ A S{iaip)] 

+ |(5c^"* A [-((/9'^?0-)o-a5 A Vv? + V9Va6V(?(TV9) + V {(p''ia)aab^ + Vip^CTab A m(p] 

= d[. . .] - if]"^ A [-5(v9^io-)o-a5V9 - 5v?Va5(io-v9) - {if^ i(y)(yab5^ - (f^(rab^{i(yv)] 

+|(5cu"^ A V[ip\iaaab + fTa6ia)v9] 
= d[. . .] + f]"^ A dl'p^iaaab + c^a6«fT)v9] + |5cu"^ A V [v?^ (ifra^b + CTa6^f^)v^] 

= d[. . .] - 1]«^ A 5[(v9V)Ca6] - |5CU"^ A V[(v9V)Ca6]. (A2) 
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We have recently discovered [|I4| that there are many identities hke ( p.2|) and ( |4.1|) in 
Riemann or Riemann-Cartan spaces of any dimension. 
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